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We study a quantum phase transition between fermion superfluid (SF) and band insulator (BI) 
of fermions in optical lattices. The destruction of the band insulator is driven by the energy gain 
in promoting fermions from valance band to various conducting bands to form Cooper pairs. We 
show that the transition must take place in rather shallow lattice height, V /Er between 2.23 and 
4.14. The latter is the prediction of mean field theory while the former is the value for opening a 
band gap. As one moves across resonance to the molecule side, the SF-BI transition evolves into the 
SF-Mott insulator transition of bosonic molecules. We shall also present the global phase diagram 
for SF-Insulator transition for the BCS-BEC family. 



There have been strong interests in ultra-cold Fermi 
gases near Feshbach resonance. Much of these interests 
have to do with the non-perturbative nature of the prob- 
lem and the richness of the phenomena. On one hand, 
this systems has universal properties at resonance[lJ 
difficult to calculate by perturbative means. On the 
other hand, the ground state turns out to be a robust 
fermion superfluid connecting the usual BCS state to the 
Bose-Einstein condensate (BEC) on different sides of the 
resonance 0. The robustness of this superfluid allows 
one to study intriguing phenomena such as those related 
to pseudo-gap and unequal spin population. Typically, 
the strength of a superfluid is reflected in its resistance 
to perturbation. The upper critical magnetic field of a 
superconductor is a good example. In the case of neu- 
tral fermion superfluids, this corresponds to the critical 
rotational frequency at which pairing is switched off. Per- 
turbation can also take the form of an optical lattice, in 
which case, the strength of the superfluid is reflected in 
the critical lattice height V^* which destroys superfluid- 
ity. Such experiment has been performed recently by the 
MIT group 

Since a homogenous Fermi gas can be changed contin- 
uously from a BCS to a BEC superfluid as a function 
of scattering length, the question is how the entire BCS- 
BEC family is affected by a lattice potential. Clearly, 
the answer depends on the fermion density. In the re- 
cent MIT experiment [3], there are two fermions per site. 
Increasing the lattice height Vb on fermion side will lead 
to a transition from BCS supefluid to a band insulator, 
whereas on the molecular side, this transition will be from 
a boson superfluid to a Mott insulator with one boson per 
site[3|. While it is the same transition, the properties of 
the superfluid and insulating phases are changed contin- 
uously across the resonance. At resonance, due to the 
universal behavior of the system, one expects the critical 
lattice height Vq (measured in recoil energy Er) at which 
superfluid-Insulator (SF-I) takes place to be a universal 
number, which is a measure of the strength of this robust 
superfluid. 

It is interesting to note that in most solid state mate- 
rial the band gap (~ leV) is several orders of magnitude 



larger than the typical pairing energy (~ 0.1 — lmeV). 
The situation is very different in strongly interacting 
Fermi gases. First of all, the pairing energy is compa- 
rable to Fermi energy. Secondly, the band gap is highly 
tunable. This offers one a unique opportunity to study 
the competition between pairing and lattice effects. 

In this letter, we shall study the transition between su- 
perfluid and insulator near Feshbach resonance by vary- 
ing lattice height. We shall determine the critical lattice 
height Vq /Eh and the size of the insulating region in a 
trap; and discuss the global phase diagram at the end. 
Our study is based on mean field theory. Since the effects 
of pair fluctuation are weak in the BCS regime but be- 
come increasingly strong in the BEC regime, mean field 
theory is only accurate on the BCS side of the resonance, 
less accurate near resonance, and becomes very poor on 
the BEC side. The situation is identical to calculating 
T c for the crossover family, where fluctuation corrections 
grow rapidly on the BEC sideQ. The simplicity of mean 
field field, however, makes it a useful first step to study 
the phase diagram on the BCS side, provided one re- 
alizes its limitations. On the other hand, we shall see 
that the mean field theory (in the presence of lattices) 
provides a surprisingly good result at resonance as it is 
sufficiently close to a rigorous lower bound. Moreover, 
the qualitative feature provided by the mean field theory 
is informative enough for one to construct a global phase 
diagram that extends to the BEC side. 

(A) Band Structure: The single particle Hamiltonian 
of a 3D cubic optical lattice is TL = hi, i = 1,2,3, 
hi = -h 2 df /(2M) + V sm 2 (Kx l ), K = tt/A, and A is the 
lattice spacing fixed by the wave length of the laser, and 
M is the mass of the fermion. Typically, Vq is measured 
in unit of "photon recoil energy" £7r = h 2 K 2 /(2M). It is 
numerically straightforward to calculate the Bloch func- 
tion and energy of the ID lattice hamiltonian h, which 
we denote as <f> n ,k and e n> fe, where n is the band in- 
dex and k is the crystal momentum. Bloch's theorem 
implies the expansion (j)n,k(x) — J2 m u n ^{m)e l ^ k+mKS>x , 
where k lies in the first Brillouin Zone (BZ). The Bloch 
state of the 3D lattice H is nk = J2 G '"nk(G)e i ( k + G ) r , 
where keBZ, nk (x) = nLi &nk, {%%), u n k(G) = 
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Ili=i u niki (nii), n = {nx,n,2,n^} is the band index, and 
G = {miK,ni2K,m3K} is the reciprocal lattice vector. 
The energy eigenvalues are e n k = J2i=i e rn,k z - 

Note that except for the (lowest) s-band (0,0,0), all 
higher bands are degenerate. For example, the p-bands 
(1,0,0), (0,1,0) and (0,0,1) are three-fold degenerate. 
Numerically, we found that only when Vq/Er > 2.23, 
a band gap E gap will develop between the s-band and 
the p-bands for this cubic lattice, as shown in FigQjA). 
FigOjB) shows how the band gap E gap grows with lattice 
height Vq. 

(B) Mean- field Hamiltonian. The many-body Hamil- 
tonian is H = H Q + Hint, H = J^nka^^^n^a, 
where V'nka crea tes fermion with quantum numbers n, k 
and the spin <r, £ nk = e n k — /i, /i is chemical potential, 
and 

Wint = ff5Z a "^m.k+PT^m' -k+p^n'.-k'+piV'n.k'+pT, 

n 

(1) 

where r\ denotes the set {n, n', m, m', k, k', p}. As the 
ultraviolet behavior of particle interaction will not be af- 
fected by lattice potentias with typical wavelength A, the 
renormlization scheme of the coupling constant g remains 
to be 1/g = m/ (4irTi 2 a s ) — l/fi E q m/h 2 q 2 , where a s is 
the s-wave scattering length, and f2 is the volume. Here, 
the q-sum is over all wavevector q, not only in the first 
BZ. 

Pairing with zero crystal momentum is described by 
terms with p = in Eq. (TT]). The corresponding coeffi- 
cient a n is then 



t ( r )^'-k( r )^',-k'(r)^ n ,k'(r) 
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where <2„ n , k = Eq u n,k( G - Q)u n / _ k (G + Q), and Q 
is a reciprocal lattice vector. The reduced BCS Hamilto- 
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FIG. 1: (A) Density of states (DOS) for V = 2E R (dashing 
line, no band gap) and Vo ~ 3Er (solid line, finite band gap) 
of a cubic lattice. -E gap is the band gap between the lowest 
band and the first excited band. (B) Band gap as a function 
of lattice depth. 



nian then assumes the form 
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construction, H m -f reduces to the non-interacting Hamil- 
tonian Ho as all Ag vanish. The insulator described by 
this mean field theory is therefore taken by assumption 
as a band insulator. This is the limitation of the theory 
and will have to be corrected in a more accurate theory. 

Within this mean field theory, the superfluid to insu- 
lator (SF-I) transition is a consequence of the competi- 
tion between the band insulator and a BCS superfluid. 
Starting from a band insulator, it is necessary for the 
pairing interaction to be strong enough to overcome the 
band gap, so that a pair of fermions from the lowest band 
can be promoted to a higher band to facilitate number 
fluctuation necessary for pairing in different bands. Al- 
ternatively, as the band gap increases, it will be more ad- 
vantageous for a superfluid with density appropriate to a 
set of filled bands to eventually collapse into these filled 
bands to minimize its energy through the band gaps. 

( C) SF to Insulator transition. Since we are interested 
in the onset of pairing order, we can expand the free- 
energy in powers of Ag, which gives 
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where Wg is the pairing susceptibility. At zero- 
temperature, it is 



w G = q- 1 E E <^ |2[0( ^ k)+e( «"'- k)1 

keBZ n,n' 



£n,k+£n', -I, 
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where ©(a;) = l/2(— 1/2) for x > 0(x < 0). The system 
is SF if at least one of Ag is not zero, which requires at 
least one Wg > —m/{A-Kh 2 a s ). 

Since u„_ k (-Q) = < k (Q), we write Q^ n , k = J2q 
<<,k(-Q-G) u n> k(G-Q)= Jd 3 r <^ k (r) n , k (r) 
e~ l2Gr . We also find numerically that Wg=o is always 
larger than Wg^q ■ So the Cooper instability first occurs 
at the zero-momentum channel. Since Q^n^k = $n,n'j 
the condition for superfluidity can be greatly simplified 
as 



Wg ~° n ( £ 2 |£ nk | £ n 2 k 2 /m ) 
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Eq.© shows that Wg=o must be of the form 

W G =o = K s E^W{n/E R , Vo/Er), (6) 
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where W is a dimensionless function. The condition 
Eq.([5]) for superfluidity then becomes 

W(n/E R , V /E R ) > -l/(8irKa s ). (7) 

The function W has the following properties: 

(a) It diverges whenever /i lies in a band, and is finite 
when fi is within a band gap. This is a consequence 
of the integrand |Cn,k| _1 = l e n,k — A* | 1 in Eq.®. The 
behavior of W for different V is plotted as a function /x 
in Fig [21 where we have only displayed the part of W for 
\x within the first band gap[6|. 

For given Vo , as /j, rises from the top of the s-band to 
the bottom of the p-band, W drops from a divergent value 
to a finite minimum and then rises back to a divergent 
value. This minimum will occur at fi*(Vo) somewhere 
inside the band gap. This minimum value of W will be 
denoted as W*(V ) = W(n*(V )/E R ,V /E R ), which is 
important for determining the phase diagram. 

(b) As V D increases, the band gap is widened. At the 
same time, as shown in Figj2j the curve W as well as its 
minimum /x*(Vo) are shifted to higher of Vo, while the 
value of W*(Vb) drops. For given a s , there is an "onset" 
lattice height (denoted as V^nset) where the minimum 
of W barely touches -l/(8irKa a ), i.e. VV*(V on sct) = 
W(fi*(V onsct )/E R ,V /E R ) = -l/(8nKa,). The chemi- 
cal potential where this touching occurs will be denote 
as /ionsct = M*(K)nset). The word "onset" means that this 
is the minimal value of Vo for band insulator to occur in 
some range of chemical potential. 

(D) Determination of phase boundary: The phase 
boundary on the fi — Vo plane separating the superfluid 
and the band insulator is given by W(fi/E R , Vo/E R ) = 



W V =4Er V = 4.5E r V = 6E r 




FIG. 2: The function W denned in Eq.© as a function of 
chemical potential I-i/Er, for Vo = 4Er, 4.5-Er and 6-Br, 
respectively. If W is above the horizontal red dashing line 
whose value is —l/(8irKa s ), the system is a superfluid for all 
chemical potentials. Insulator occurs when a portion of W 
falls below — l/(87rifo s ). In the case of Vo = 6Er, this occurs 
in the chemical potential interval 
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FIG. 3: (A) The fi — Vo phase diagram for different scattering 
length, where —l/(Ka 3 ) = 1, and —1 respectively. The 
dashing line denotes the bottom of the lowest band, which 
monotonically increases as the increase of zero-point energy. 
The dotted vertical line T is the trajectory of local chemical 
potential form the center of the trap (c) to the edge of the 
cloud (e). (B) Solid line: Vonsct as a function of — l/(Ka s ). 
Dashing line: the lowest lattice depth for opening up a band 
gap. (C) Schematic density profile in trap for V > Vonsct and 
V < Vonsct respectively. 



— l/(8nKa s ), which can be determined graphically as 
shown in Fig. 2. As discussed before, we shall focus on 
l/(Ka s ) < 1, i.e. on the BCS side and around resonance 
but not on the BEC side. The graphical method works 
as follows: 

(i) If W intersects —l/(8irKa s ) so that it lies below 
— 1/(8ttK a s ) over an interval of /Lt, say, (ne, fi u ), as shown 
in the case of Vo = 6£"r in Fig. 2, the system is a band 
insulator within this interval of /i. See also Fig[3] 

(ii) As Vb decreases, the curve W rises. As a result, 
the interval (fi(.,fj, u ) shrinks, and the values fig and \x u 
get closer to each other. When Vo reaches a critical 
value Vo nso t, the insulating interval shrinks to zero, and 
fit = Mit = Monsct, as shown in the case of Vo = 4.5-Er . In 
this way, the points jit and [i u trace out a phase boundary 
shown in FigJ3] These values Vo nsc t and Vonsct are identi- 
cal to those defined in point (b) in Section (C). In Fig. 3, 
the phase boundaries for different l/(Ka s ) have been dis- 
placed. The onset lattice height Vonsct as a function of 
a s has also shown in FiglH^B). At resonance, it shows 
that Vonset = 4.14i?R. Since a band gap only appears for 
Vo > 2.23-Er, the mean field theory therefore predicts 
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FIG. 4: Schematic of how the fi — Vo phase diagram evolves 
across a Feshbach resonance, from the molecule side (up one) 
to resonance regime (middle one) to the fermion side (bottom 
one). 



that in a three dimensional cubic lattice, the critical lat- 
tice height for the onset of the first insulating phase is 
between 2.23-Er and 4.14_Efj at resonance. It should be 
noted that these numbers predicted here depend on the 
lattice type (such as f.c.c. and b.c.c), since they have dif- 
ferent band structure. It should also be noted that the 
lattice height A.UEr for SF-I transition is rather shal- 
low. At this lattice height, tight binding approximation 
is not valid and also the system can not be described by 
a single band attractive Hubbard model, 
(iii) In the presence of a trap, within the local density 
approximation, the chemical potential is a function of 
position, fi(r) = fi — V(r), where V(r) is the trapping 
potential. Going from the center to the surface of the 
gas, a vertically downward trajectory T is generated in 
the fi-Vo space. If this trajectory passes through the in- 
sulating region as shown in Fig|3]JA), this will show up 
as a step structure in the in-situ density profile as shown 
in FigHC). 

Starting from a superfluid, as Vq increases, an insulat- 
ing lobe will move toward the physical trajectory T from 
the right in Figj3^A). When the lobe touches T, which 
happens when Vq = V onsc t with the point of contact oc- 
curs at Vonsct, the insulating phase will develop at the 
location r where /x(r) = /i nset- 

(E) Global phase diagram: It is useful to consider the 



phase diagam in the space of chemical potential fx, lattice 
height Vo, and interaction parameter —l/(Ka s ). In the 
absence of lattice, the ground state can evolve continu- 
ously from a BCS state to a BEC by shrinking the size of 
the pair wave function continuously Q- In the case of two 
fermion per site, as V increases, the fermion superfluid 
on the BCS side (which are made up of Cooper pairs of 
Bloch states) will evolve continuously into a bosonic su- 
perfluid with one boson per site on the average on the 
BEC side. For sufficiently high V , the ground state on 
the BCS side will be a band insulator, whereas it will be 
a Mott insulator (or a Fock state) with one boson per site 
on the BEC side. Since a band insulator can be written 
as a Fock state with two fermions in the same Wannier 
state on each site. A band insulator can also evolve con- 
tinuously into a Mott insulator by simple changing the 
wavefunction of a fermion pair with opposite spin on the 
same site. In this way, there is a crossover of both su- 
perfluid and insulator as one crosses the resonance. Note 
that while filling the lowest (s) band on the BCS side 
requires two fermions per site, filling an additional (p) 
band requires totally eight fermions per site. The Mott 
insulators on the BEC side, however, can have any in- 
teger number of bosons, or any even integer number of 
fermions per site. To go from the BEC side to the BCS 
side, only Mott phases with 1 and 4 bosons can crossover 
to a band insulator. Mott phases with 2 and 3 bosons (or 
4 and 6 fermions) per site must disappear as it approaches 
the fermion side. The phase diagram will therefore look 
like FigH] schematically. 
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